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Abstract. The rich acoustic oscillation spectrum in solar-type variables make these stars particularly interesting for studying 
fluid-dynamical aspects of the stellar interior. I present a summary of the properties of solar-like oscillations, how they are 
excited and damped and discuss some of the recent progress in using asteroseismic diagnostic techniques for analysing low- 
degree acoustic modes. Also the effects of stellar-cycle variations in low-mass main-sequence stars are addressed. 
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INTRODUCTION 

Solar-type stars possess extended surface convection 
zones. The observed oscillation modes generally behave 
as acoustic modes and their frequencies are sensitive pre- 
dominantly to the sound speed c in the stellar interior. It 
appears that all possible oscillation modes are intrinsi- 
cally stable. They are excited stochastically by the strong 
emission of acoustic noise by the turbulent velocity field 
in the upper convectively unstable layers of the star. The 
excitation occurs in a broad frequency range, giving rise 
to a rich pulsation spectrum. The amplitudes of the os- 
cillations are small, typically 5 ppm L /M (Kjeldsen & 
Bedding 1995). Such small amplitudes allow us to de- 
scribe the pulsations with linear theory, however, they 
also pose a challenge to the detection limits of observ- 
ing campaigns. But the recently made progress in con- 
tinuously lowering the detection threshold of today's in- 
strumentation has allowed us to have access to these low- 
amplitude pulsators, providing high-precision data such 
as from the space mission CoRoT (Convection Rotation 
and planetary Transits; Baglin 2003), from the recently 
launched Kepler satellite (Basri et al. 2005, Christensen- 
Dalsgaard 2007), and from the ground-based observing 
campaigns, such as SONG (Grundahl et al. 2007) and 
those of the kind organized by Kjeldsen, Bedding and 
their colleagues (e.g. Kjeldsen et al. 2005). Solar-like 
oscillations have been convincingly observed on other 
stars and augur asteroseismic diagnosis that will raise 
stellar physics to a new level of sophistication. It will 
not be possible to measure either the internal structure 
or the internal motion of a distant star with a resolution 
comparable with that which we have achieved for the 
Sun, because high-degree modes will not be accessible 
in the foreseeable future; but some theoretically impor- 
tant properties, such as the gross structure of the energy- 
generating core and the extent to which it is convective, 
and the large-scale variation of the angular velocity, will 
become available. Such information will be of crucial 



importance for checking, and then calibrating, the the- 
ory of the structure and evolution of stars, the backbone 
of theoretical astrophysics. 
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FIGURE 1. Hertzsprung-Russell diagram of stars in which 
solar-like oscillations have been detected (adapted from Aerts 
et al. 2008; the original figure was kindly provided by Fabien 
Carrier). 

The first indication of excess acoustic power in other 
stars with a frequency dependence similar to the Sun was 
reported by Brown et al. (1991) for the F5 star Procyon 
A (a CMi), in which the first unambiguous detection 
of solar- like oscillations was reported by Martic et al. 
(1999), confirmed later by Mosser et al. (2008). First (un- 
confirmed) detection of individual peaks in the acoustic 
power spectrum from high-precision time-resolved spec- 
troscopic observations was published for the G5 star r] 
Boo by Kjeldsen et al. (1995), but it was not before 
2003 that an unambiguous confirmation was established 
by Carrier et al. (2003) and Kjeldsen et al. (2003). The 
location of stars in the Hertzsprung-Russell diagram, in 
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FIGURE 2. Small section of a solar acoustic power spectrum. The radial order n and spherical degree / are indicated in pairs of 
(n,l) for each mode. The large and small frequency separations, Av n j and 5v n j are in general functions of n and I and can be used 
to infer the mass and age of a star (adapted from Christensen-Dalsgaard 2001). 



which solar-like oscillations have been detected until to- 
day, are indicated in Fig.Q] and a summary of such stars 
was presented recently by Bedding & Kjeldsen (2007). 

SOLAR-LIKE OSCILLATION 
PROPERTIES 

Only modes of low degree can be observed in distant 
stars, however, in general, the observed modes are of 
high radial order, which allows us to extract diagnos- 
tic properties of the frequencies v„ / with radial order n 
and spherical degree / for the asymptotic limit n/l — > °°. 
The diagnostic properties of this type of mode have 
been studied extensively in the solar case. From asymp- 
totic theory we find for the cyclic oscillation frequen- 
cies (Gough 1986, 1993; see also Tassoul 1980 and Van- 
dakurov 1967) 
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is the inverse of twice the sound travel time between the 
centre and surface (R is surface radius), and 
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The frequency-dependent coefficient e is determined by 
the reflection properties of the surface layers, as is the 
small correction term B, and L 2 = 1(1 + 1). The value 



of Vo can be estimated from taking the average (over n 
and /) of the so-called large frequency separation Av„ / = 
V n l — between modes of like degree and consecu- 
tive order. The last two terms on the right-hand side of 
equation (Q~|) lift the degeneracy between modes with the 
same value of n + 1/2 and leads to the so-called small fre- 
quency separation 8v n j = V„ / — V n -\,i+2- One obtains a 
frequency structure in which modes of odd degree fall 
approximately halfway between modes of even degree, 
which is illustrated in Fig. [2] for a solar spectrum. The 
mean small frequency separation (averaged over n) 



(5v„,,) = (v„,-v„_ M+2 ) ~ [2A(2/ + 3)]- 
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is predominantly determined by the acoustic sound speed 
in the stellar core and hence is sensitive to the chemical 
composition there and consequently is an indicator for 
the stellar age (e.g., Gough 2001, Houdek & Gough 
2008). 

Seismic diagnostic 

Acoustic modes depend predominantly on the sound 
speed in the stellar interior and consequently on the 
chemical composition. As the star evolves, the chemical 
composition will vary and consequently also the value of 
the small frequency separation 8v n j. So does the mean 
large frequency separation (averaged over n and I) 

(Av) = (v„, - v„_ u ) ~ vo - (M/fl 3 ) 1 / 2 ( 5 ) 

mainly because of the increasing surface radius with 
time. A convenient way to demonstrate the dependence 
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FIGURE 3. Stellar evolutionary tracks (solid curves, indi- 
cated by the stellar mass) and curves of constant central hy- 
drogen abundance (dashed curves, indicated by the central hy- 
drogen abundance by mass, X c ) in terms of the average large 
separation (Av) and average small separation (Sv) = {8v„q) 
[c.f. equations $5$ & ©] (from Christensen-Dalsgaard 1993). 



of the sound speed and radius on the effects of stellar 
evolution can be illustrated in a ((Av),(5v)) diagram, 
depicted in Fig. [3] 

Additional information of structural aspects of solar- 
like stars is obtained from the seismic signatures con- 
tained in Av„./ and 5v„./. Abrupt variation in the stratifi- 
cation of a star (relative to the scale of the inverse radial 
wavenumber of a seismic mode of oscillation), such as 
that resulting from the (smooth, albeit acoustically rela- 
tively abrupt) depression in the first adiabatic exponent 
7= (d]np/d]np) s caused by the ionization of helium, 
where p, p and s are pressure, density and specific en- 
tropy, or from the sharp transition from radiative to con- 
vective heat transport at the base of the convection zone, 
induces small-amplitude oscillatory components (with 
respect to frequency) in the spacing of the cyclic eigen- 
frequencies v„ / of seismic oscillation and consequently 
also in Av„./ and 8v„j. We call such abrupt variations an 
acoustic glitch. One might hope that the variation of the 
sound speed c induced by helium ionization might enable 
one to determine from the low-degree eigenfrequencies 
a measure that is directly related to, perhaps even almost 
proportional to, the helium abundance, with little con- 
tamination from other properties of the structure of the 
star. 

A convenient and easily evaluated measure of the os- 
cillatory component produced by acoustic glitches is the 
second multiplet-frequency difference with respect to or- 
der n amongst modes of like degree /: 

A 2 v„,/ = v„_i,/-2v„., + v„ +M (6) 
(Gough 1990). Any localized region of rapid variation 
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FIGURE 4. The symbols in the upper panel show second 
differences A2V n j (see equation [6j for low-degree (/ = 0, ..,3) 
modes obtained from simulated data for a 1 M© model of age 
5.54 Gy. The simulation was based on two periods of 4-months 
observations, separated by 1 year, with the SONG network (see 
Grundahl et al. 2007). The solid curve is a fit to A 2 v based on 
the analysis by Houdek & Gough (2007). The dashed curve is 
a smooth contribution, modelled as a third-order polynomial 
in V , which represents near-surface effects. The lower panel 
displays the remaining individual contributions from the acous- 
tic glitches to A2V: the dotted and solid curves show the con- 
tribution from the first and second stages of helium ionization 
and the dot-dashed curve is the contribution from the acoustic 
glitch at the base of the convective envelope. 



of either the sound speed c or the density scale height, 
or a spatial derivative of them, induces an oscillatory 
component in A 2 v (from here on the subscripts n,l have 
been dropped for simplicity) with a 'cyclic frequency' 
approximately equal to twice the acoustic depth 

,R 

X= c _1 dr (7) 

■'''glitch 

of the glitch, and with an amplitude which depends on 
the amplitude of the glitch and which decays with v once 
the inverse radial wavenumber of the mode becomes 
comparable with or less than the radial extent of the 
glitch. 

Various approximate formulae for the oscillatory com- 
ponents that are associated with the helium ionization 



have been suggested and used, by e.g., Basu et al. (1994, 
2004), Monteiro & Thompson (1998, 2005) and Gough 
(2002), not all of which are derived directly from explicit 
acoustic glitches. Gough used an analytic function for 
modelling the dip in the first adiabatic exponent. In con- 
trast, Monteiro & Thompson assumed a triangular form. 
Basu et al. have adopted a seismic signature for helium 
ionization that is similar to that arising from a single 
discontinuity; the artificial discontinuities in the sound 
speed and its derivatives that this and the triangular repre- 
sentations possess cause the amplitude of the oscillatory 
signal to decay with frequency too gradually, although 
that deficiency may not be immediately noticeable within 
the limited frequency range in which adequate astero- 
seismic data are or will imminently be available. More 
recently Houdek & Gough (2007) proposed a seismic di- 
agnostic in which the variation of y in the helium ion- 
ization zone is represented with a pair of Gaussian func- 
tions. This correctly results in a decay of the amplitude 
of the seismic signature with oscillation frequency that 
is faster than that which the triangular and the single- 
discontinuity approximations imply, and also takes some 
account of the two ionization states of helium. Moreover, 
they incorporated the acoustic cutoff frequency into the 
variation of the eigenfunction phase with acoustic depth, 
thereby improving the discrepancy between the seismi- 
cally inferred depths of the acoustic glitches and that of 
a corresponding stellar model. An example of the appli- 
cation of Houdek & Gough's technique is presented in 
Fig.|U which shows the resulting fit to simulated data for 
a solar-like star. 

Mode parameters 

Were solar p modes to be genuinely linear and sta- 
ble, their power spectrum could be described in terms 
of an ensemble of intrinsically damped, stochastically 
driven, simple-harmonic oscillators, provided that the 
background equilibrium state of the star is indepen- 
dent of time (Fig. [5]); if we assume further that mode 
phase fluctuations contribute negligibly to the width of 
the spectral lines, the intrinsic damping rates of the 
modes, Y/2, could then be determined observationally 
from measurements of the pulsation linewidths F. 

The power (spectral density), P, of the surface dis- 
placement %ni(t) of a damped, stochastically driven, 
simple-harmonic oscillator, satisfying 
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which represents the pulsation mode of order n and de- 
gree /, with linewidth T„/ and frequency co„i and mode 
inertia I B \, satisfies 
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FIGURE 5. Power spectral density of a randomly excited, 
damped, harmonic oscillator. Pf represents the spectral density 
of the random force and P is the product of the Lorentzian P^ 
and Pf (adapted from Kosovichev 1995). 



assuming T,,/ <C (0„i, where /(f) describes the stochastic 
forcing function. Integrating equation (O over frequency 
leads to the total mean energy in the mode 
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where A„i is the displacement amplitude (angular brack- 
ets, (), denote an expectation value) and V„j is the rms 
velocity of the displacement. The total mean energy of a 
mode is therefore directly proportional to the rate of work 
(also called energy supply rate) of the stochastic forcing 
Pf at the frequency a),,/ and indirectly proportional to F n [. 

Equations Cf8T> — (TT~0t are discussed in terms of the dis- 
placement 4 (from here on we omit the subscripts n and 
/), but in order to have a direct relation between the ob- 
served velocity signal v{t) = d<i; /dt and the modelled ex- 
citation rate Pf we shall first take the Fourier transform 
V(v) of v(t) (v = co/2n). It follows that the total mean 



energy E of the harmonic signal of a single pulsation 
mode is then given by (Chaplin et al. 2005) 



E = IV 2 = 18 
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in which 
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is the maximum power density - which corresponds to 
the 'height' of the resonant peak in the frequency do- 
main (see Fig.|2]i. The height H is the maximum of the 
discrete power, i.e. the integral of power spectral density 
over a frequency bin 5 = 1 /T ^ s , where r o t, s is the total 
observing time. The following expressions 
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and 



H := 
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(172)2/ 

provide a direct relation between the observed height H 
(in cm 2 s 2 Hz '), the modelled energy supply rate Pf (in 
erg s~'), and damping rate rj = Y/2. 

PULSATION COMPUTATIONS AND 
AMPLITUDE RATIOS 

The linearized pulsation equations for nonadiabatic ra- 
dial oscillations can be presented as (e.g. Balmforth 



1992a): 
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where 8 is the Lagrangian perturbation operator, and 
for simplicity the right hand side of the perturbed mo- 
mentum equation is formally expressed by the func- 
tion / (the full set of equations can be found in, e.g., 
Balmforth 1992a). Equations $15[ are solved subject to 
boundary conditions to obtain the eigenfunctions and 
the complex angular eigenfrequency (0 = (Or + irj, where 
(Or is the (real) pulsation frequency and rj =Y/2 is the 
damping rate in (s _1 ). The turbulent flux perturbations 
of heat and momentum, 8L C and 8pu and the fluctu- 
ating anisotropy factor 54> are obtained from the non- 
local, time-dependent convection formulation by Gough 
(1977a,b). 

From the linearized nonadiabatic pulsation equations 
( fT3T > theoretical intensity-velocity amplitude ratios 
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TABLE 1. Absorption lines and their wavelengths X of vari- 
ous helioseismic instruments. Also listed are the optical depths 
T5000 at 5000 A and the corresponding approximate heights 
above the photosphere (h = at T = T e ff) at which the lines 
are formed. 



Instrument 


line 


A (A) 


^5000 


height (km) 


BBSO 


Ca 


6439 


0.05 


~ 129^ 


BiSON 


K 


7699 


0.013 


~25oJ 


MDI 


Nil 


6708 


9xl0~ 3 


^300^. 


GOLF 


Na D1/D2 


5690 


5xl0~ 4 


~ 500" 



* from Libbrecht (1988) 

+ from Christen-Dalsgaard & Gough (1982) 

** from Toutain et al. (1997, but see also Baudin et al. 2005) 



can be compared with observations, without the need of 
a specific excitation model and all its uncertainties in 
describing the turbulent spectra. 

It is, however, important to realize that various instru- 
ments observe in different absorption lines and conse- 
quently at different heights in the atmosphere. This prop- 
erty has to be taken into account not only when com- 
paring observations between various instruments (e.g. 
Christensen-Dalsgaard & Gough 1982), but also when 
comparing theoretical amplitude estimates with observa- 
tions (Houdek et al. 1995). Table 1 lists some of the rel- 
evant properties of various instruments. 

In the top panel of Fig. [6] the theoretical amplitude 
ratios (equation ( fT6] l) of a solar model are plotted as a 
function of height for several radial pulsation modes. The 
mode energy density (which is proportional to rp l l 2 8r) 
increases rather slowly with height; the density p, how- 
ever, decreases very rapidly and consequently the dis- 
placement eigenfunction 8r increases with height. This 
leads to the results shown in the upper panel of Fig. [6] 
where the decrease in the amplitude ratios with height is 
particularly pronounced for high-order modes for which 
the eigenfunctions vary rapidly in the evanescent outer 
layers of the atmosphere. It is for that reason why so- 
lar velocity amplitudes from, e.g., the GOLF instrument 
have larger values than the measurements from the Bi- 
SON instrument (by about 25%, Kjeldsen et al. 2005). 
The lower panel of Fig. [6] compares the estimated so- 
lar amplitude ratios (curves) with observed ratios (sym- 
bols) as a function of frequency. The model results are 
depicted for velocity amplitudes computed at different 
atmospheric levels. The observations are obtained from 
accurate irradiance measurements from the IPHIR in- 
strument of the PHOBOS 2 spacecraft with contempora- 
neous low-degree velocity data from the BiSON instru- 
ment at Tenerife (Schrijver et al. 1991). The thick solid 
curve represents a running-mean average, with a width of 
300 /iHz, of the observational data. The theoretical ratios 
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FIGURE 6. Top: Calculated amplitude ratios (see equation 
J16I Q as a function of height in a solar model for modes with 
different frequency values. Bottom: Theoretical amplitude ra- 
tios (surface luminosity perturbation over velocity) for a solar 
model compared with observations by Schrijver etal. (1991). 
Computed results are depicted at different heights above the 
photosphere (fr=0km at T = T e ff) The thick, solid curve indi- 
cates a running-mean average of the data (from Houdek et al. 
1995). 



for h = 200 km (dashed curve) show reasonable agree- 
ment with the observations. 

In Fig. [7] model results for the F5 star Procyon A are 
compared with observations (horizontal dotted line) by 
Arentoft et al. (2008). Theoretical results are shown for 
two stellar atmospheres: a VAL-C (Vernazza et al. 1981) 
atmosphere scaled with the model's effective tempera- 
ture T e ff (thin curves), and for an Eddington atmosphere 
(thick curves). For both stellar atmospheres the agree- 
ment with the observations is less satisfactory than in the 
solar case, indicating that we may not represent correctly 
the shape of the pulsation eigenfunctions. Consequently 
there is need for adopting more realistically computed 
atmospheres in the equilibrium models, particularly for 
stars with much higher surface temperatures than the 
Sun. It should, however, be mentioned that the current 
photometric observations are still uncertain. 
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FIGURE 7. Calculated amplitude ratios (see equation d!6l >) 
for a model of Procyon A are compared with observations by 
Arentoft et al. (2008; horizontal dotted line). Theoretical results 
are shown for a scaled VAL-C atmosphere (black, thin curves) 
and for an Eddington atmosphere (red, thick curves). Top: The 
theoretical amplitude ratios are shown as a function of height 
and for three different pulsation modes. The frequencies of the 
three pulsation modes are indicated. Bottom: The theoretical 
amplitude ratios are shown as a function of frequency at three 
different heights in the stellar atmosphere. The heights above 
the photosphere h = km are indicated. 



DAMPING RATES 

Damping of stellar oscillations arises basically from two 
sources: processes influencing the momentum balance, 
and processes influencing the thermal energy equation. 
Each of these contributions can be divided further ac- 
cording to their physical origin, which was discussed in 
detail by Houdek et al. (1999). 

Important processes that influence the thermal energy 
balance are nonadiabatic processes attributed to the mod- 
ulation of the convective heat flux by the pulsation. This 
contribution is related to the way that convection mod- 
ulates large-scale temperature perturbations induced by 
the pulsations which, together with the conventional K- 
mechanism, influences pulsational stability. 

Current models suggest that an important contribution 
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FIGURE 8. The symbols in the left-hand panel are the measured linewidths A,,/ = T„i /2n (we denote the FWHM in unit of 
cyclic frequency by A„j) of solar low-degree p modes obtained from a 3456-d data set collected by BiSON between 1991 and 2000 
(Chaplin et al. 2005). The data are compared with the theoretical damping rates jt^'t],,/ (connected by the solid curve) obtained 
from the model computations discussed in the previous section (from Chaplin et al. 2005). In the right-hand panel theoretical results 
of (2k) ~ 1 r\ n [ (solid curve) by Dupret et al. (2004) are compared with observations of A,,/ / 2 (symbols). 



that influences the momentum balance is the exchange 
of energy between the pulsation and the turbulent ve- 
locity field through dynamical effects of the fluctuating 
Reynolds stress. In fact, it is the modulation of the turbu- 
lent fluxes by the pulsations that seems to be the predom- 
inant mechanism responsible for the driving and damp- 
ing of solar-type acoustic modes. It was first reported 
by Gough (1980) that the dynamical effects arising from 
the turbulent momentum flux perturbation Sp t contribute 
significantly to the damping F. Detailed analyses (Balm- 
forth 1992a) reveal how damping is controlled largely 
by the phase difference between the momentum pertur- 
bation and the density perturbation. Therefore, turbulent 
pressure fluctuations must not be neglected in stability 
analyses of solar-type p modes. 

A comparison between the latest linewidth measure- 
ments (full- width at half-maximum) A„i = T„i/2n and 
theoretical damping rates is given in Fig. [8] The observa- 
tional time series from BiSON (Chaplin et al. 2005) was 
obtained from a 3456-d data set and the linewidths of the 
temporal power spectrum extend over many frequency 
bins 8 = l/2r o b s . In that case the linewidth in units of 
cyclic frequency is related to the damping rate according 
to 

A„i = 7t- l ri„, . (17) 

Recently Dupret et al. (2004) performed similar sta- 
bility computations for the Sun using the time-dependent 
mixing-length formulation by Gabriel et al. (1975, 1998) 
and Grigahcene et al. (2005), which is based on the for- 
mulation by Unno (1967). Their results are illustrated in 
the right panel of Fig. [SJ which also shows the character- 
istic plateau near 2.8 mHz. It is, however, interesting to 
note that their findings suggest the fluctuating convective 



heat flux to be the main contribution to mode damping, 
whereas for the model results shown in the left panel of 
Fig. [8] which are based on Gough's (1977a,b) convection 
formulation, it is predominantly the fluctuating Reynolds 
stress that makes all modes stable. 

Houdek et al. (1999) computed damping rates r\ of 
solar-like oscillations in about 160 stars with masses be- 
tween 0.9 Mq and 2.0 in the vicinity of the main 
sequence, and for various metallicities and convection 
parameters. Recently Chaplin et al. (2009) conducted a 
similar theoretical study of estimating mode lifetimes 
T = TJ in low-mass stars and compared the theoreti- 
cal estimates with the latest linewidth measurements in 
twelve solar-type variables. They found that the mean 
mode lifetimes of the five most prominent solar-like p 
modes (t) scale like 

(T)-r eff 4 . (18) 

Their results are depicted in Fig. [9] where the diamond 
symbols are the theoretical estimates of (t) obtained 
from a grid of models computed in the manner of Chaplin 
et al. (2005) and the crosses with error bars are the 
linewidth measurements of twelve main-sequence, sub- 
giant and red-giant stars. 

STOCHASTIC EXCITATION 

Because of the lack of a complete model for convection, 
the mixing-length formalism still represents the main 
method for computing the turbulent fluxes in the convec- 
tively unstable layers in a star. One of the assumptions in 
the mixing-length formulation is the Boussinesq approx- 
imation, which results in neglecting the acoustic wave 
generation by assuming the fluid to be incompressible. 
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FIGURE 9. Average mode lifetimes (t) (averaged over the 
five most prominent p modes) in low-mass stars. The diamond 
symbols are the results from the stellar equilibrium and pulsa- 
tion calculations and the cross symbols are observations of 12 
main-sequence, sub-giant and red-giant stars. The solid curve 
represents the power law , where T e g is the surface temper- 
ature of the models that were used in the stability computations 
(from Chaplin et al. 2009). 



Consequently a separate model is needed to estimate the 
rate of the acoustic noise (energy supply rate) generated 
by the turbulence. The excitation process can be regarded 
as multipole acoustic radiation (Lighthill 1952). Acous- 
tic radiation by turbulent multipole sources in the context 
of stellar aerodynamics has been considered by Unno 
& Kato (1962), Moore & Spiegel (1964), Unno (1964), 
Stein (1967), Goldreich & Keeley (1977), Bohn (1984), 
Osaki (1990), Goldreich & Kumar (1990), Balmforth 
(1992b), Goldreich, Murray & Kumar (1994), Musielak 
et al. (1994), Samadi & Goupil (2001) and Chaplin et al. 
(2005). 

The mean amplitude A of a mode is determined by a 
balance by the energy supply rate Pf from the turbulent 
velocity field and the thermal and mechanical dissipation 
rate characterized by the damping coefficient 77 [see 
equation (I10H . The procedure that we adopt to estimate 
A is that of Chaplin et al. (2005), whose prescription 
follows that of Balmforth (1992b). 

We represent the linearized pulsation dynamics by the 
simplified equation 



dt 2 



■-F(u) + G(s') (19) 



for the displacement | (r, t ), which is now also a function 
of radius r, of a forced oscillation corresponding to a 
single radial mode satisfying the homogeneous equation 



(20) 



in which co (and £) are real and JC is a linear spatial 
operator. The (inhomogeneous) fluctuating terms on the 



right-hand-side of equation dT9l > arise from the fluctuat- 
ing Reynolds stresses 



F(u) = V • (puu - (puu)) 



(21) 



and from the fluctuating gas pressure (due to the fluctuat- 
ing buoyancy force), represented by G(s'), where s' is the 
Eulerian entropy fluctuation (Bohn 1984; Osaki 1990; 
Goldreich & Kumar 1990; Balmforth 1992b; Goldreich, 
Murray & Kumar 1994; Samadi & Goupil 2001). The 
latest numerical simulations by Stein et al. (2004) sug- 
gest that both forcing terms in equation ( fT9l ) contribute 
to the energy supply rate Pf by about the same amount, 
a result that was also reported by Samadi et al. (2003) 
using the turbulent velocity field and anisotropy factors 
from numerical simulations (Stein & Nordlund 2001). In 
this paper we consider only the term of the fluctuating 
Reynolds stresses and because we use only radial modes, 
only the vertical component F3 of F is important, 



F i (u i )~—(pul-{pul)). 
or 



(22) 



If we define the vertical component of the velocity cor- 
relation as P33 = (M3M3), its Fourier transform P33 can 
be expressed in the Boussinesq-quasi-normal approxi- 
mation (e.g. Batchelor 1953) as a function of the turbu- 
lent energy spectrum function E(k, a>): 



P33 — 



link 2 



(23) 



where £ is a wavenumber and is an anisotropy param- 
eters given by 

¥ = [24>/3(4>-l)] 1/2 , (24) 

which is unity for isotropic turbulence (Chaplin et al. 
2005). This factor was neglected in previously published 
excitation models but it has to be included in a consistent 
computation of the acoustic energy supply rate. Follow- 
ing Stein (1967) we factorize the energy spectrum func- 
tion into E(k,co) = E(k)Q.((o;Ti ( )^ where x k = X/kuk is 
the correlation time-scale of eddies of size %jk and ve- 
locity Uk, the correlation factor X is of order unity and 
accounts for uncertainties in defining Tfc. 

The energy supply rate is then given by (see Chaplin 
et al. 2005 for details) 

n r R ( dB \ 2 



with 



y(r,(o) 



fc- 2 £ 2 (K-)n(T fc ;o))dfc, (26) 



where K = kt/n, t is the mixing length, R is surface ra- 
dius, and ^ r is the normalized radial part of % . The spec- 
tral function 5? accounts for contributions to Pf from the 
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FIGURE 10. Top: Predicted apparent velocity amplitudes 
(defined to be v2 times the rms value) for a model of a Cen A, 
computed according to equation < | 1 3 1 > . Bottom: Linear damping 
rates for a model of a Cen A, obtained by solving the fully 
nonadiabatic pulsation equations l |15| >. 



FIGURE 11. Top: Predicted apparent velocity amplitudes 
(defined to be y2 times the rms value) for a model of £ Hydrae, 
computed according to equation J 1 3b . Bottom: Linear damping 
rates for a model of t, Hydrae, obtained by solving the fully 
nonadiabatic pulsation equations J 1 5 1 > (adapted from Houdek 
& Gough 2002). 



small-scale turbulence and includes the normalized spa- 
tial turbulent energy spectrum E{k) and the frequency- 
dependent factor Q.(T k ;co). For E(k) it has been com- 
mon to adopt either the Kolmogorov (Kolmogorov 1941) 
or the Spiegel spectrum (Spiegel 1962). The frequency- 
dependent factor £l(Xk', (O) is still modelled in a very rudi- 
mentary way and we adopt two forms: 
- the Gaussian factor (Stein 1967) , 



£2 G ((o;Tt) 



2k 



(27) 



- the Lorentzian factor (Gough 1977b; Samadi et al. 
2003; Chaplin etal. 2005), 



1 



7rV21n2 1 + (coz k /y/21n2) 



(28) 



The Lorentzian frequency factor is a result predicted for 
the largest, most-energetic eddies by the time-dependent 
mixing-length formulation of Gough (1977b). Recently, 
Samadi et al. (2003) reported that Stein & Nordlund's 



hydrodynamical simulations also suggest a Lorentzian 
frequency factor, which decays more slowly with depth z 
and frequency co than the Gaussian factor. Consequently 
a substantial fraction to the integrand of equation d25l ) 
arises from eddies situated in the deeper layers of the 
Sun, resulting in a larger acoustic excitation rate Pf. 

OSCILLATION AMPLITUDES 

Model predictions of the mode height H were computed 
according to equation ([Pil l. As in the solar case damp- 
ing rates for other stars are obtained from solving the 
eigenvalue problem ( fTBI ) and the energy supply rates are 
calculated from expression (f25j). With these estimates 
for 77 and Pf Houdek & Gough (2002) predicted ve- 
locity amplitudes for several stars, using equation dT3T >. 
Results for stochastically excited oscillation amplitudes 
and linear damping rates in the solar-like star a Cen A 
and in the sub-giant % Hydrae are illustrated in FigsfTUl 
andQT] Kjeldsen et al. (2005) reported mode lifetimes for 
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FIGURE 12. Predicted apparent velocity amplitudes (de- 
fined to be \fl times the rms value) as function of light-to-mass 
ratio for stochastically excited oscillations in other stars. Obser- 
vations from several authors are plotted by the plus and triangle 
symbols. The theoretical estimates by Houdek & Gough (2002) 
are plotted as diamond symbols. The scaling law of Kjeldsen 
& Bedding (1995) is illustrated by the dotted line and results 
reported by Samadi et al. (2005) are indicated by the dashed 
line. 



a Cen A of about 2.1 days (but see also Fletcher et al. 
2006), which are in reasonable agreement with the the- 
oretical estimates of about 1.7 days for the most promi- 
nent modes (the mode lifetime z= rp 1 ; see lower panel 
of Fig. [Toll. For ^ Hydrae, however, the theoretical mode 
lifetime of the most prominent modes is about 18 days 
which is in stark contrast to the measured value of about 
2-3 days by Stello et al. (2004, 2006), yet the estimated 
velocity amplitudes for E, Hydrae are in almost per- 
fect agreement with the observations by Frandsen et al. 
(2002). 

A comparison between predicted and observed veloc- 
ity amplitudes in several solar-type pulsators is illus- 
trated in Fig. [T2] The dotted line is the scaling law by 
Kjeldsen & Bedding (1995), which is based on the com- 
putations by Christensen-Dalsgaard & Frandsen (1983), 
and the dashed line is the scaling relation reported by 
Samadi et al. (2005) using the convective velocity pro- 
files from 3D numerical simulations (Stein & Nordlund 
2001), a Lorentzian frequency factor in equation ( 126) , 
and the theoretical damping rates from Houdek et al. 
(1999). For the cooler stars the theoretical results (scaling 
laws) are in reasonable agreement with the observations, 
whereas for hotter stars, such as for Procyon A, the theo- 
retical velocity amplitudes are overestimated by both the 
scaling laws and the stochastic excitation models. 

Recently, Chaplin et al. (2009) suggested a new scal- 
ing law for the intensity amplitudes by combining their 
finding of the mean mode lifetime, (t) °^ T e f^ [cf. equa- 
tion with Kjeldsen & Bedding's (1995) scaling law 



for the intensity amplitudes inferred from a narrow-band 
observation of wavelength A, i.e. {8L/L)x °= {L/M) /T^, 
leading to 

H « g" 2 (29) 

for the maximum mode height. Since the surface gravity 
changes fairly slowly along the main sequence, the new 
scaling relation (|29l for H, which assumes narrow-band 
intensity observations, suggests that stars notably cooler 
than the Sun might have mode heights that are compa- 
rable to those solar-like pulsators that are hotter than the 
Sun. Moreover, Hekker et al. (2009) reported that the in- 
tensity heights in about 780 red giant stars, observed in 
broadband photometry with the CoRoT satellite, follow 
the scaling law H °< g~ 2 2 , a result that supports the the- 
oretical finding d29]l by Chaplin et al. 2009. 

STELLAR-CYCLE EFFECTS 

From helioseismic data, such as those provided by the 
BiSON, we have learnt that not only the oscillation fre- 
quencies change with time over the 11 -year solar cycle 
but also the height H and width Y. In the solar case 
the absolute fractional change from solar activity min- 
imum to maximum in (H) (angular brackets () denote 
an average over the five most prominent p modes) is 
about 40%. That in (r) is about 20% (Chaplin et al. 
2000). We therefore expect also in solar-like oscillators 
not only the mode frequencies but also the mode heights 
and linewidths to show stellar-cycle variations. 

Chaplin et al. (2007, 2008) estimated stellar-cycle fre- 
quency shifts and mode height variations for a grid of 
31 models with masses between 0.7 M and 1.3 M and 
stellar ages, f*, in the range from the ZAMS to 9Gy, 
using the Padova isochrones (e.g., Bonatto et al. 2004) 
to specify mass, radius, effective temperature, r e ff, and 
luminosity. The evolutionary tracks of these models are 
presented in the left panel of Fig. Qj] in terms of a (?*, 
Teff) plot. To predict the stellar-cycle changes Chaplin 
et al. used the Ca II H&K index for surface activity in 
stars. This index is usually expressed as R' HK , the aver- 
age fraction of the star's total luminosity that is emit- 
ted in the H&K line cores. The authors then used the 
data of 22 main-sequence stars collected by the Mount 
Wilson Ca II H&K programme (e.g., Saar & Branden- 
burg 2002) from which the R' HK cycle amplitude val- 
ues, A/?y K , were determined. The A/?^ K values were then 
simply scaled against the 0.4 jj,Hz frequency shift seen 
for the most prominent low-/ modes in the Sun, assum- 
ing that the frequency shifts 5v cyc scale approximately 
linearly with AR' HK . In order to estimate the frequency 
shifts of the 3 1 grid models, it was necessary to calculate 
first the R' HK values for the 31 models. This was done ac- 
cording to the procedure by Noyes (1983) and Noyes et 
al. (1984). The R' HK values so obtained were then used 
to estimate the corresponding cycle frequency shifts by 
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FIGURE 13. Left: Age - effective temperature plot for a sequence of stellar models, obtained from the Padova isochrones, for 
which the stellar-cycle induced frequency shifts and mode height variations were predicted. Right: Predicted stellar-cycle frequency 
shifts (left ordinate) and relative height variations (right ordinate) for a sequence of stellar models with masses between 0.7 M s and 
1 .3 M (from Chaplin et al. 2008). 



interpolating linearly in the (<5v cyc , ^hk) table of the 22 
main-sequence stars observed by the Mount Wilson Ca II 
H&K programme (for a detailed discussion see Chaplin 
et al. 2007). The outcome of this procedure for estimat- 
ing the stellar-cycle frequency shifts of the most promi- 
nent p modes in the stellar models is depicted in the right 
panel (left ordinate) of Fig. [13] Also shown in this figure 
are the predicted fractional changes in (H) (right ordi- 
nate), obtained in a similar way as the frequency shifts 
and assuming a fractional change of (H) of 40% for the 
Sun. These results suggest that, in the considered age 
range of 1 - 9 Gy, the variation of p-mode frequencies 
and heights over the stellar cycle can be up to 1.5-2 
times larger than in the Sun, and that these variations de- 
pend predominantly on stellar age and less on r ett - or stel- 
lar mass. 
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